The anisotropic glassy properties of quasicrystals within an amended tunnehng model 
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We use an extended version of the standard tunneling model to explain tlie anisotropic sound 
absorption in decagonal quasicrystals. The glassy properties are determined by an ensemble of two 
level systems (TLS), arbitrarily oriented. The TLS is characterized by a 3 x 3 symmetric tensor, 
which couples to the strain field through a3x3x3x3 tensor of coupling constants, [[i?]]. The 
structure of [[i?]] reflects the symmetry of the quasicrystal. 

We also analyze the probability distributions of the elements of [T] in this particular model for a 
better understanding of the characteristics of "isotropic" and "anisotropic" orientations distributions 
of the ensemble of TLSs. We observe that the distribution of the elements is neither simple, nor 
intuitive and therefore it is difficult to guess it a priory, using qualitative arguments based on the 
symmetry properties. 

PACS numbers: 



INTRODUCTION 



Despite of almost four decades of study of the glassy 
properties materials [iHll , the nature of the two- level sys- 
tems (TLS), the ubiquitous hypothetical microscopic en- 
tities that are held responsible for these properties, is still 
not known. As a general picture, it is accepted that they 
represent dynamic defects, which are atoms or groups 
of atoms that tunnel from one minimum energy config- 
uration to another. But these atoms are in general not 
identified and, even in the cases when they are identified 
(like in crystals with defects), the TLS spectrum cannot 
be obtained based on a microscopic model. 

The study of solids with anisotropic glassy properties 
is especially interesting since this brings additional infor- 
mation about the TLSs and requires a critical perspective 
on the standard tunneling model (STM). 

To identify the origin of the TLS the thermal and 
acoustical properties of a number of crystalline systems 
with defects |4l-ll2l| and quasicrystals jl3l4l5l | at low tem- 
peratures have been investigated thoroughly. As ex- 
pected, the glass like properties in both the disordered 
crystalline system and in the quasicrystals were revealed. 
However, a crucial difference in comparison with ordinary 
amorphous solids which possess TLS excitations was de- 
tected. Namely, a pronounced anisotropy in internal fric- 
tion was clearly marked. As a consequence of that, the 
baffling physical problem about the possible origin of the 
anisotropy appeared. Indeed, there are two competitive 
ways for the explanation. First, the effect is explained by 
an anisotropy in the distribution over the "orientations" 
of the TLSs in the ensemble. The interaction between the 



TLS and the elastic field is described by the Hamiltonian 



where 5 = 2jijSij and [S] is the strain field of the 
phonon-we assume everywhere summation over the re- 
peated indices. The symmetric second rank tensor [7] 
characterizes the TLS and its "deformability" under elas- 
tic strain. The anisotropy of the physical properties is a 
refiection of the values taken by the elements of [7] , which 
are determined by the lattice symmetries. Bert, Bcl- 
lessa and Grushko made a conjecture regarding these val- 
ues from which they eventually recovered the anisotropic 
sound attenuation rates in decagonal quasicrystals |15| . 

In the second approach [IGMlSj -the one that we shall 
employ in this paper-the TLS is characterized by a 3 x 3 
symmetric tensor [T] and the coupling between [T] and 
[S] is made through a forth rank tensor of coupling con- 
stants denoted by [[R]]. Explicitly, jij = RkUjTki and 
5 = [T]* : [[i?]] : [S]. The elements of [T] are determined 
by a unit vector t, which is the direction of the TLS, 
wheres the structure of [[i?]] is determined by the sym- 
metries of the host material. In this model, even if the 
TLSs are isotropically oriented, the anisotropy of the sys- 
tem is imposed by the properties of the tensor of coupling 
constants. 

We applied the second approach to the crystalline 
materials of different symmetries with embedded TLSs, 
assuming the isotropy of the orientations of the TLSs 
H-iij and we calculated the attenuation of ultrasound 
waves of different polarizations and propagating in differ- 
ent crystallographic directions. Not only that this model 
describes in a simple way the asymmetries of the glassy 
properties, but it also allows us to make predictions about 
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the relative attenuation rates of sound propagating in dif- 
ferent directions. 

In this paper we apply the model of Ref. [l^ to the at- 
tenuation of ultrasound waves in quasicrystals and we ob- 
tain the attenuation rates along different crystallographic 
directions expressed in terms of the components of the 
tensor of coupling constants. We a pply our calculations 
to the experimental results of Ref. [15| , which enable us 
to calculate the relations between some of the compo- 
nents of the coupling constants tensor. 

In order to better understand the distribution of the 
TLSs in isotropic and anisotropic materials, we calculate 
the distribution of the elements of [T] for an isotropic 
distribution of unit vectors t. This distribution is neither 
simple, nor trivial, and by this we draw attention to the 
fact that it is difficult to guess it a priori, without a more 
detailed model. 



THE DESCRIPTION OF THE TLS 

Let us introduce the notations by presenting briefly the 
model. The Hamiltonian of the free TLS is 
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where A is called the asymmetry of the potential and A 
the tunnel splitting. The eigenvalues of i?TLS B^re ±e/2, 
where e = a/A^ -f A^ is the excitation energy of this TLS. 
The groundstatc will be denoted by 1 1) and the excited 
state by | t). The interaction Hamiltonian of the TLS 
with the strain field is given by ([IJ 

The parameters A and A are distributed with the prob- 
ability P(A,A) = Pq/^, where Pq is a constant. If ex- 
pressed in terms of e and it = A/e, the probability distri- 
bution becomes P{e, u) = Po/{u\/l ~ u^). 

As usual, we work in the abbreviated subscript no- 
tations and write [S] and [7] as six-dimensional vec- 
tors: S = (511,^22, ^33, 2^23, 2^13, 25i2)* and 7 = 
(711)722,733,723,713,712)*, where the superscript t de- 
notes the transpose of a matrix or a vector. 

As stated before, jij = RkUjTki- Assuming that the 
TLS is characterized by a direction in space, t, the tensor 
T is formed of the components of t [16[ . In abbreviated 
subscript notations T = (t^, t^., 2tytz,2tztx, 2txty)'^ , 
Rijki becomes and 7 = [i?]-^ • T [16[. The structure 
of [R] is determined by the symmetries of the lattice, 
since 5 is a scalar and should be invariant under coor- 
dinates transformations, whereas [R] should be invariant 
under the symmetry transformations that leaves the lat- 
tice invariant [l6| . 

The absorption rate of a phonon-with wavenumber k 
and polarization cr-by a TLS is 



rk<.(t) = ^ A!^|T* . [R] ■ Sk.P<5(e - huj). (3) 



The main characteristic of the TLS-phonon interaction is 
contained in the quantity Mk,cr(t) = T* • [R] ■ Sktr, which 
bears an intrinsic anisotropy through the matrix [R], on 
which the symmetries of the lattice are imposed. 

The average scattering rate of a phonon by the ensem- 
ble of TLSs is obtained by averaging rkcr(t) over e. A, 
and t. In this way we get the total phonon absorption 
rate, 

^k.' = Y^"(l^'^k.(t)P>nk.tanh (^^^ , (4) 
which may be put into the standard form, 

= Y^oll^^'kW. tanh , (5) 



where N = yjfij (2V piS) is the normalization constant of 
the phonon mode and 7,^^^ = (|Mka(t)n^/^/(7VA:) is the 
(average) phonon-TLS coupling constant. 

We have no model regarding the orientations of the 
TLSs, so we shall assume that t is isotropically oriented. 

The structure of \R\ should be similar to that of the 
elastic stiffness constants, [c], so for a decagonal qua- 
sicrystal should have the form 15, 2^ 
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where rge = {r\\ — ri2)/2-the axis z is taken along the 
tenfold axis. We observe that the structure of both, \R\ 
and [c] , is similar in decagonal quasicrystals and in hexag- 
onal lattices. 

Like in hexagonal lattices, in the decagonal quasicrys- 
tals we can have pure longitudinal and transversal waves 
propagating in all the three directions, a;, y, and z. 

The coupling constants, 7^. ^ are similar to those cal- 
culated for hexagonal lattices in pH : 
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(7a) 
(7b) 
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where by I and t we refer to longitudinal and transversal 
polarizations, respectively. 

While for the longitudinal waves the direction of polar- 
ization is obvious, for the transversal waves the direction 
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of polarization is indicated by the second unit vector in 
the subscript of 7, in the Eqs. (|7c)) and ((7d]) . 

Due to the isotropy condition in the decagonal plane, 
the coupling constants of the phonons propagating in this 
plane are independent of the direction of propagation if 
they have similar polarization. 

Compariiig now Eqs. ((7c)) and (j7dp with the results of 
Bert et al. uM, we obtain 



ni - ri2 
2r44 



4.2, (8) 



where and are obvious notations from Ref. [l^l ■ 
To find the ranges of the coupling constants, we rewrite 
the Eqs. (I7a|) and ([Tbl as 
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(9a) 
(9b) 



Since the function f{a,b) = + 5^ + 2a6/3 satisfies 
/(o, b) > /(0, 0) = for any a and 6, we find from 
Eqs. dH) that 



llij^lly,i>^ and 7L,/> " 
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(10) 



The minimum values of the coupling constants, 7^^ ^ and 
7I2 for given ri3 (fTO)) . are obtained for rn = ri2 = 
— ri3/4 and = — 2r33/3. The ratio between these 
minimum values is 

(7L.i)mi„/(7fey,i)mi„ = 3/8. (11) 

and is independent of rn, ri2, ri3, and r33. 



case in the STM. To understand the difficulties to define 
by general, qualitative arguments the "anisotropy" of 7 
in the STM, we calculate in this section the probability 
distribution of the elements of our tensor T, under the 
assumption that the unit vectors t are isotropically ori- 
ented. We shall see that these elements do not have all 
the same distribution of probability, nor the distributions 
are constant for the ranges of these variables. 

We define t by the Euler angles, 9 and 0: 
tx = sin cos ^, ty = sin 6* sin and = cos 6*. Keeping 
tz and (j> as variables, we write T* = [{l — tD cos^ 0, (1 — 
tl) SIB^ (|3,tl,2tz^/U^Sm(|},2tz^/l^ coscj}, (1 
t2)sin(20)]. 

The variables tz and are uniformly distributed in the 
intervals [0,1] and [0, 27r], respectively. We denote their 
probability distribution by Pi{tz,4>) = l/(47r). Then the 



probability distribution in the variables Tj (/ = 1, 
and d is 
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The simplest of all is eventually 

which, if integrated over (p, gives 
^^3(^3^ ^ 
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For the components Ti and T2 we obtain the probabil- 
ity distributions 
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THE ANISOTROPY OF THE GLASSY 
PROPERTIES 

The anisotropics of the glassy properties of a disor- 
dered system are determined by the properties of the 
tensor 7. Bert, Bellessa and Grushko made in Ref. [l^ 
a conjecture regarding the values of 7,^ , which, they con- 
sidered, would lead to the observed anisotropy. 

The distributions we obtained in this paper are differ- 
ent from those conjectured in Ref. [l3|- In our model 
the "orientations" of the TLSs are well defined by the 
unit vectors t and the symmetries of the host material 
are incorporated into the coupling constant tensor, [R\. 
Therefore the anisotropic distribution of the elements of 
the tensors T and 7 are clearly defined. This is not the 



where 0{x) is the Heaviside step function. An extra factor 
of 2 appears in the Eqs. because Ti and T2 are even 
functions of tz- 

For an isotropic distribution of t, the components Ti, 
T2, and arc equivalent and using Eqs. (fT4|) and (fT5|) 
we obtain the identity 
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(16) 

The components T4, T5, and Tq are also equivalent for 
isotropic t. For example 
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FIG. 1: The probability distributions of the components of 
the tensor T: Pi{T) = P2{T) = P:i{T) = l/(2^/T), for T G 
(0,1] (solid line), and 2P4(|r|) = 2Ps,{\T\) = 2P,i{\T\), for 
T e [-1, 0) U (0, 1] (dashed line). 



the tensor T - which describes the TLS - under the as- 
sumption that the directions of the ensemble of TLSs, 
defined by the unit vectors t, arc isotropically oriented. 
We observe that the distributions of the elements of T 
are neither (all six) similar to each other, nor they are 
constant distributions. In conclusion, this distribution 
is not possible to find only by qualitative, general argu- 
ments, but one needs for this a deeper model, like the 
oneofRef. 
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where r6/sin(20) = 1 - > o. From Eq. P7a|) we get 



P6(T) = 1- f (17b) 
- P4(T) = P5(r), (17c) 

where T G [-1, 1] and P4(|T|) = Psd^l) = PedPl). 

We notice that limT-j-i PxeiT) — l/(2-\/2), whereas in 
0; PTe{T) has a logarithmic divergence. 

The equalities (|17cP are obtained using the equivalence 
between the components T4, T<^, and Tg. 

We observe that the isotropic distribution of the TLS 
orientations do not correspond to a constant distribution 
of the values of the components of T, nor to equal distri- 
butions these values. Therefore it is not straightforward 
to draw conclusions regarding the relations between the 
distributions of the values of 7/, based only on general 
arguments. 



CONCLUSIONS 

In this paper we describe the anisotropy of the glassy 
properties of the decagonal quasicrystals in the the model 
of Ref. • We show that the glassy properties of these 
quasicrystals are similar to those of hexagonal disordered 
lattices l2l| and we obtained the TLS-phonon average 
coupling constants, 7^^., which arc dependent on the 
phonon's propagation direction, k, and polarization, a. 
We apply the results to the experimental data of Ref. |15| 
and we obtain the ratio {ree/r 44)'^ « 4.2, where 7-44 and 
rge are components of the tensor of coupling constants. 

In order to better understand the characteristics of 
"isotropic" and "anisotropic" distributions of TLSs, we 
calculate the probability distributions of the elements of 
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